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Abstract. We study efficient combinatorial algoritlims to produce tlie 
Hasse diagram of tiie poset of bounded faces of an unbounded polyfie- 
dron, given vertex-facet incidences. We also discuss the special case of 
simple polyhedra and present computational results. 



1. Introduction 

The bounded subcomplex of an (unbounded) convex polyhedron, is its set 
of bounded (or, equivalently, compact) faces, partially ordered by inclusion. 
Polytopal complexes of this type arise in several situations each of which is 
interesting in its own right. Prominent examples are the tight spans of finite 
metric spaces of Dress see also Isbell |16| . and the tropical polytopes of 
Develin and Sturmfels [9] . The purpose of this note is to present algorithms 
to deal with such objects. 

Typically, in applications a bounded complex of a polyhedron P is given 
only implicitly by a set of inequalities defining P. The primary goal is to 
establish algorithms to make the bounded subcomplex explicit from this 
input. A direct approach is to enumerate the full face lattice £j{P) of a 
polytope P projectively equivalent to P (which exists if P is pointed) and 
to filter it to obtain B{P), the poset of bounded faces. However, since 'C(P) 
often is much larger than B{P), this is not efficient. 

One natural approach is to start with a (dual) convex hull computation 
which also yields the vertex-facet incidences. Here we focus on combina- 
torial algorithms which take these vertex-facet incidences as input. Other 
approaches are discussed briefly. Our Algorithm [1] in Section 13.11 is a mod- 
ification of a combinatorial algorithm for face lattice enumeration |20j . It 
uses the vertex-facet incidences of P. This method of selective generation is 
best possible, in the sense that its running time is linear in the size of B{P). 
If only the vertex-facet incidences of P are given, that is, no information 
about the unbounded edges is present, one can still compute B{P), although 
at the cost of a higher running time (quadratic in the size of B{P)). The 
corresponding Algorithm [2] is based on interleaving the computation of the 
bounded faces with an incremental computation of the poset's Mobius func- 
tion. 

For the sake of simplicity of exposition we usually consider complexity 
questions in the RAM model. However, it is easy to modify each of our 
results such that the bit complexity can be determined: The only potential 
source of non-polynomiality in terms of RAM complexity arises from calling 
LP type oracles. Therefore, whenever necessary, we explicitly mention the 
relevant sizes of the linear programming problems which need to be solved. 
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If additional structural information is available, specialized algorithms 
come into focus. For instance, in the case where the polyhedron P is simple, 
it turns out that B{P) is determined by the (directed) vertex-edge-graph 
of P, see |15| Section 2]. This is an unbounded version of a result of Blind 
and Mani [5] obtained by applying techniques due to Kalai |21| . In Section U 
we employ the reverse-search approach of Avis and Pukuda [3] to generate 
the (directed) vertex-edge graph of B{P) from an inequality description of P, 
provided that P is simple. This can be used to either construct the bounded 
faces or to efficiently compute the face numbers of P and B{P). 

In Section [5l computational experiments with Algorithm [1] are presented. 
We investigate five different cases in which it is interesting to compute the 
bounded subcomplex. The computations show that some (surprisingly) large 
instances can be handled. The paper closes with a list of open problems. 

We are grateful to an anonymous referee for meticulous reading and for 
requiring a correction in Algorithm [2j This also contributed to a cleaner 
description of this method. 



2. Preliminaries and Notation 

Let P be a polyhedron. The lattice of all faces {face lattice) of P is denoted 
by £(P). We define ip{P) = |£(P)|, that is, the number of all faces of P. 
Moreover, let n be the number of vertices of P and m be its number of facets. 

A polyhedron is called pointed if it does not contain an entire affine line. 
In the sequel we always assume that P is pointed but unbounded. For 
our purposes this does not mean any loss of generality, since the bounded 
subcomplex of a polyhedron is non-empty if and only if it is pointed. In this 
case, the polyhedron P is projectively equivalent to a polytope P. Each such 
polytope is called a projective closure of P. A projective closure of P has a 
special face Poo (the far face) corresponding to the face at infinity. Fixing 
an admissible projective transformation 7 which maps P to P, the face Pqo 
is the unique maximal face among the faces of P that are not images of 
faces of P under 7. Note that the combinatorial type of Pqo depends on the 
geometry of P, not only its combinatorial structure, and that its dimension 
can be any number in {0, . . . , dimP — 1}. 

We begin with the description of an algorithm to compute the polytope P 
from an inequality description of P C W^. We refer to Ziegler's mono- 
graph |26| for a general discussion of admissible projective transformations 
and |19l §3.4] for an explicit construction. First, we compute one vertex v 
of P by solving a linear program similar to Phase I of the Simplex Algo- 
rithm. Such a vertex exists as P is pointed. In a second step we determine 
the active constraints at v, that is, those inequalities which are satisfied with 
equality at v. Let r be the affine transformation moving v into the origin. 
In the image t(P) the constraints active at correspond to homogeneous 
linear equations. Any subset of those constraints defines a polyhedral cone 
containing the translated polyhedron t(P). Among the constraints active 
at we choose a dual basis by Gau£- Jordan-elimination. Next we pick a 
linear transformation p mapping this dual basis to the basis which is dual 
to the standard basis of W^. Then the image p[t{P)) is a polyhedron with 
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vertex such that the boundary hyperplanes of the positive orthant define 
vahd inequahties. In particular, p{t{P)) is contained in the positive orthant. 

In general, intersecting the image of a polyhedron in R'^ under a projec- 
tive transformation in PGLrf(R) with M'^ may yield something non-convex. 
However, if M G GLrf_|_i(M) is an invertible matrix with non- negative coef- 
ficients, then the induced projective linear transformation [AI] € PGL(i(R) 
maps the positive orthant into itself, and the image of any polyhedron inside 
is again a polyhedron. Now let fj, be the unique projective linear transfor- 
mation which fixes each coordinate hyperplane of M*^ and which additionally 
maps the hyperplane at infinity to the affine hyperplane ^ = 1 . This 
transformation can be written as a (d -|- 1) x (d -|- l)-matrix of zeroes and 
ones. The polyhedron 

P := Mt{P))) 

is contained in the simplex given by ^ Xj < 1 and Xi > for i £ [d]. In 
particular, P is bounded and projectively equivalent to P. The equation 
^ Xj = 1 defines a hyperplane supporting P, and its intersection with P is 
the far face Foq. For given matrix A and right hand side b, let i{A, b) be the 
time complexity of solving the linear feasibility problem to find a point x 
satisfying Ax < b. The discussion above can be summarized as follows. 

Proposition 1. If P = P{A,b) is given by inequalities, there exists an 
explicit algorithm to compute P in 0{£{A,b)) time. 

The unbounded edges of P correspond to edges of P that contain exactly 
one vertex in F^o . The vertices of P correspond to vertices of P that are not 
contained in Foo- Via these relations it is easy to describe the face lattice 
of P in terms of P and vice versa. The following relationship holds between 
the size ip = (p{P) of the face posets of P and the size Ip = (p{P) of its 
projective closure. 

Lemma 2. We have Tp < 2{ip — 1). 

Proof. For each unbounded face F oi P let s{F) be the intersection of F 
with the hyperplane at infinity. Then s is a map from the set of unbounded 
faces of P to the set of non-empty faces of the far face F^o which is surjective. 
This proves that P has at most twice as many non-empty faces as P. To 
establish our slightly stronger claim, observe that P must have at least one 
vertex, since we assumed P to be pointed. □ 

The above bound is tight: For P consider any pyramid with its basis as 
its special face. In view of the obvious inequality < ^ the lemma shows 
that Tp G ©(</?); hence in statements about asymptotic complexity, <p and Tp 
can be used interchangeably. 

We label the facets of P from 1 to m and identify each facet with its 
index in F := {1, . . . ,rn}. Similarly, label the vertices of P from 1 to n, and 
identify each vertex with its index in V := {1, . . . Let / G {0, l}"^^" be 
a vertex-facet incidence matrix of P with entries I{f,v), that is, I{f,v) = 1 
if facet / contains vertex v, and /(/, v) = otherwise. Denote by a the 
number of incidences, that is, the number of ones in /. The face lattice of 
a polytope is atomic and co-atomic: Every face of P is uniquely defined by 
its set of vertices and the set of facets it is contained in, respectively. In the 
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following, we will usually identify a face with its set of vertices and store each 
such set by its characteristic vector (bitset). Hence, the storage requirement 
per face is in 0(n). Modifications for other formats are straightforward - 
see also Section [5] below. 

The computation of an incidence matrix for a polytope given by inequali- 
ties requires to generate the set of vertices, which can be exponentially many 
in the number of inequalities. There are several algorithms to generate the 
vertices of a polytope, see Seidel |23| for an overview and Avis et al. [2] as 
well as |17| for discussions of the complexity of vertex generation algorithms. 
In fact, it is currently unknown whether the vertices of a polytope can be 
generated in polynomial time in the combined size of the input and out- 
put. It is, however, known that it is NP-complete to decide whether the list 
of vertices of a polyhedron is complete (unbounded edges are ignored), see 
Khachiyan et al. |22j and Boros et al. [7]- For simple d-polytopes, in which 
each vertex is contained in exactly d facets, the reverse search algorithm by 
Avis and Fukuda [3] produces the vertices in time 0{d -m -n). 

Let T> be the Hasse diagram of P. That is, T> is a directed rooted acyclic 
graph whose nodes correspond to the elements of If Nh, Nq are nodes 

in T> and H and G are the corresponding faces of P, then there is an arc 
{Nh,Ng) in V if and only if H C G and dim(G) = dim(i?) + 1. Note 
that we distinguish between a face F and its corresponding node Np in D, 
because this makes a difference in Algorithms [1] and [2] below. 

3. The Face Poset of the Bounded Subcomplex 

The partially ordered set (poset) of the bounded faces of an unbounded 
polyhedron P is denoted as B{P), and we call it the bounded subcomplex of P. 
As a polytopal complex, B{P) is not pure in general: it may have maximal 
faces of various dimensions. The bounded subcomplex is non-empty if and 
only if P is pointed. We denote the size of B{P) by ip' in the following. 

We primarily address the problem to compute B{P), where P or P is 
given in terms of a vertex-facet incidence matrix. Other input is discussed 
at the end of this section. The output B{P) should be given by its Hasse 
diagram, which is just a subgraph of the Hasse diagram of P. Additionally 
or alternatively, the nodes, that is, faces, might be labeled by either the 
set of facets the corresponding face is contained in or by its set of vertices 
(this representation is unique, since the faces are bounded); this requires an 
additional overhead of at least 0{n) or 0[m) per face, respectively. 

For the case where ip' k, Lp and a vertex-facet incidence matrix / S 
{0, 1}"*^" of P is available, the direct approach is to apply the algorithm 
of |20| to generate the Hasse diagram of P in 0(n • a • (^) time if the faces 
are represented by their vertices; this algorithm will be modified in the next 
section. We then remove the unbounded faces and their incident arcs by 
checking whether the intersection of each face with is empty or not; if 
the faces are stored as bitsets, this requires 0(n) time per face; the total 
running time for this step is then dominated by the generation of the Hasse 
diagram of P . 

Alternatively, one can also generate the Hasse diagram by working on the 
dual, using 0(m •«•(/?) time; in this case, faces are represented by the facets 
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1 
2 
3 
4 
5 
6 
7 
8 
9 
10 

11 



Input: incidence matrix I of P, far face Fqo C V 
Output: Hasse diagram P of B{P) 
initialize T) with corresponding to the empty face 
initialize the list Q C V{V) x 2^ by (iV0,0) 
while Q 7^ do 

choose some {Nh, H) G Q and remove it from Q 
compute the set Q of incident faces G with dim G = dim H +1 
foreach G G ^ do 

if G is bounded then 

locate/create the node Nq corresponding to G in D 
if Ng was newly created then 
\_ add {Ng,G) to Q 

add the arc {Nh, Nq) to V 
Algorithm 1: The method of selective generation 



they are contained in. In this case, the removal step takes 0{a) time per 
face, since we have to intersect the vertex sets of F^o and of the facets that 
contain the face (storing the vertex sets in sorted lists). If the intersection 
is empty, the face is bounded. Thus, its running time is again dominated by 
the generation of the Hasse diagram. 

li ip' <^ (p there is a more efficient algorithm available, which we describe 
in the next section. 

3.1. Selective Generation. In the following, we will describe an algorithm 
that computes (the Hasse diagram of) B{P) efficiently if an incidence matrix 
of P is available. Our Algorithm [1] is a simple modification of the algorithm 
in |20) . which computes the face lattice of a polytope. We provide some 
details for completeness and to enable a running time analysis. 

The algorithm performs a graph search through the Hasse diagram of B{P) 
starting from the bottom (the empty face). We denote by V{T>) the set of 
nodes of D. In each step we consider a face H that has not been considered 
before, see Step HI Then we generate the set Q of all faces G with G D H, 
dimG = dimH + 1, that is, there exists an arc {Nh,Ng) between the nodes 
Nh and Nq corresponding to H and G, respectively. It is shown in |20j 
that the generation of Q can be performed in O(n^) C 0(n • a) time. For 
each G G ^, we can now test whether G is unbounded by checking whether 
its intersection with F^o is empty or not (in 0{n) time), see Step [71 If it 
is bounded, we proceed to generate the arc {Nh,Ng). For this, one has to 
find or create the node Nq corresponding to G (Step [8|). Here, we use a 
data structure, called face tree in |20| . to store (bounded) faces and their 
corresponding nodes. 

Each edge of a face tree stores a vertex. The vertices along a path from 
the root form a generating set for each face. Each sub-path corresponds to 
a sub-face. Because each face of a bounded face is bounded as well, each 
node of the tree corresponds to a different bounded face. The key point is 
that the computation of the face from the generating set takes 0(a) time, 
see |20| for more details. For future reference we state the following. 
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Lemma 3 (^20j). The face tree data structure allows to store bounded faces 
such that finding a face (or determining whether the face is not present) takes 
0{a) time with a storage requirement of 0{n ■ k), if k faces are stored. 

Returning to Algorithm [U the while-loop in Step [3] is executed for each 
of the if' bounded faces, because of the boundedness test in Step [71 Hence, 
Steps [H and [5] contribute 0(n • a ■ ip') time in total. The for-loop in Step [6] is 
executed for each out-arc of a face in B{P) with respect to the Hasse diagram 
of P. Since a face can have at most n out-arcs, the for-loop is executed at 
most n- (f' times. Step [7] takes 0{n) C 0(a) time and Step [8] uses 0{a) 
time, because of Lemma [S] note that since P is pointed, each facet contains 
at least one vertex, i.e., n < a. Thus, the for-loop contributes 0(n • a ■ ip') 
time. This shows the following. 

Theorem 4. Given the vertex-facet incidences of P, the Hasse diagram 
of B{P) can he computed in time 0{n ■ a ■ ip'). 

Remark 5. By Lemma [3l the face tree uses 0(n • (p') space. Storing faces 
as bitsets, the list Q needs at most 0(n • ip') space. The output amounts to 
an additional space requirement of 0(n • ip'). Thus, we need a total amount 
of storage of 0(n • (/?'). 

Remark 6. Algorithm [T] is faster than the straight-forward algorithm de- 
scribed above \{n-p'< min{n, m} • ip. 

In Section H] below, we discuss the special case of simple polyhedra. 

3.2. Selective Generation Without Knov^ing the Face at Infinity. 

If only a vertex-facet incidence matrix / G {0, 1}"*^" of an unbounded 
(pointed) polyhedron P is known, but no information about the unbounded 
edges is available, one can still produce (the Hasse diagram of) B{P) as 
follows. 

Let y{P) := {vert(F) : F proper face of P} U {0} (where vert(F) is the 
set of vertices of F) be the poset of the vertex sets of proper faces of P. It 
can be computed from any vertex-facet incidence matrix of P, since it is the 
set of all non-empty intersections of the subsets of V defined by the rows of 
the incidence matrix (and additionally the empty set). Note that this poset 
contains the poset B{P) as a subposet, but may contain additional vertex 
sets of unbounded faces. It follows that ip' < ip" , where ip" is the size of '?{P); 
see Remark [13] below for an example which shows that the gap between ip' 
and (p" can be large. Moreover, because B{P) is a polytopal complex, we have 
the following: if vert G C vert F for an unbounded face G and a bounded 
face F, there exists a bounded face H such that vert H = vert G. 

Using results of |18| . the boundedness of faces can be decided in the fol- 
lowing way. The Mobius number of a poset element S € 5'(P) is defined 
as 



where the sum ranges over all poset elements S' € '^{P) strictly less than S. 
For a face H ^ P, we define fi{H) := fi{veTtH). For H = P, we add an 
artificial top-element 1 to J'(P) and define ft{P) ■= ^(1)- Then, the face H 




if 5 = 0, 
f^{S'), otherwise, 

S'cs 
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Input: incidence matrix I oi P 
Output: Hasse diagram V of B{P) 

initialize T> with corresponding to the empty face, M(0) = {(0, 1)} 
initiahze queue Q C V{V) x 2^ by (A^0,0) 
while Q / do 

remove first {Nh , H) in Q 
compute Mobius number fi{H) using M(ff) 
if ii{H) (H is bounded) then 
foreach vertex v not in H do 

G ^ intersection of vertex sets of facets containing H U {v} 
locate/create the node Nq corresponding to G in 2? 
if Ng was newly created then 
add {Ng,G) to Q 
M(G) ^ M{H) 

M(G) ^ M(G) U {{H, fi{H))} 
add the arc {Nh, Nq) to V 

Algorithm 2: Method using the Mobius function 



is bounded if and only if fi{H) ^ 0, see |18l Corollary 4.5]. Once all poset 
elements are known, the entire Mobius function of J'(-P), that is, the Mobius 
numbers for all poset elements, can be computed by inverting an appropriate 
(if" X <^")-matrix (the so-called (^-matrix) in 0{{ip")^) time. 

We obtain the following algorithm to compute the Hasse diagram of B{P). 
Since J'(P) is atomic, its Hasse diagram can be generated by the algorithm 
in |2U| in time 0{n ■ (3 ■ ip"), where f3 is the number of vertex-facet incidences 
of the polyhedron P (excluding unbounded information). Then the vertex 
sets corresponding to unbounded faces are removed by using the Mobius 
function. In total, we obtain an 0(max{n • /?, (9?")^} ' time algorithm. 
Since 'P{P) is also co-atomic, one can again apply the algorithm of |20| to 
the dual, which then yields an 0(max{m • (3, (v?")^} ' time algorithm. 
One can often improve these running times as follows. 

Theorem 7. Given the vertex-facet incidences of P (without information 
on unbounded edges), the Hasse diagram of B{P) can be computed in time 
0(max{n^ ■ ip' ,n ■ 13} ■ n ■ ip'). 

Proof. Algorithm [2] is a modification of Algorithm [TJ in which the Mobius 
numbers jl{H) are computed on the fly in Step [5j To this end, we use a 
set JA{H) that stores all elements of y{P) strictly below H and their cor- 
responding Mobius numbers. Then Step [5] is a straight-forward summation. 
Throughout this algorithm, the faces G and H oi P are encoded as sets of 
their vertices. 

To update M(ff) correctly, we perform a breadth-first search (BFS) over 
the Hasse diagram of y{P) from bottom to top, by organizing Q in the 
algorithm as a queue. We propagate the necessary information for the com- 
putation of Mobius functions in Steps[T2]and[T3l It follows that once a face H 
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leaves Q in Step 31 'M{H) contains all elements of J'(P) strictly below H and 
their Mobius functions. 

For the correctness of the entire algorithm it is essential that fi{H) is 
computed correctly in Step [5j This is a consequence of two facts: First, 
fi{H) is computed after all elements of J'(-P) below H have been processed; 
this is different from Algorithm [TJ Second, the data structure 'M(H) may 
also contain vertex sets of unbounded faces, but their contribution to p,{H) 
is zero, by \18\ Corollary 4.5]. Therefore, fi{H) evaluates to the Mobius 
number of H in the poset y(-P). 

If H is bounded, only its covering elements G in '?{P) are generated in 
Step [HI More precisely, G is generated from H by adding one vertex v not 
contained in H and computing the closure with respect to /; that is, G is 
the intersection of all vertex sets of facets containing H U {v}; see |20| §2.2]. 
This is the only situation how new elements can enter the queue. Hence, 
the while- loop in Step [3] is executed for each of the (f' bounded faces plus at 
most n ■ ip' faces that are unbounded. 

The above discussion implies that |M(i?)| < n-ip' , and hence Step[5]takes 
0{n-ip') time. Note that the size of M(i?) is 0{n^-ip'), as we store the vertex 
set of each face as a bitset. The sets M(G) can be organized using a face-tree. 
As observed in KP Lemma [3] can be used in this context as well, where we 
have to replace a by /3 and n by n. To add information to M(G) in Step ll3l 
we need to check whether the subset H is present in M(G), which can then 
be done in time 0{n ■ j3). To copy JA{H) in Step[T2l we need 0{in? ■ ip') time. 

All other steps are as in Algorithm [1] and their analysis is similar as in the 
proof of Theorem m The only difference is that a can be replaced by /3, that 
is, it is easy to see that by using the vertex-facet incidence matrix of P, we 
only generate sets in J'(P) . In total, we obtain an 0(max{n^ • n • /3} • n • if') 
time algorithm. □ 

Remark 8. Each face-tree data structure needs 0(n^ • ip') space, see Re- 
mark [5] In total, the algorithm requires 0{n^ • (<^')^) space for maintain- 
ing the sets M(-ff). If faces are stored as bitsets, the queue Q needs at 
most 0{'n? • if') space. The output amounts to an additional space require- 
ment of 0{n ■ if'). Thus, we need a total amount of storage of 0{rfi ■ (v?')^). 

3.3. Polyhedra Given in Terms of Inequalities. As mentioned above, 
the bounded subcomplex B{P) can be obtained by removing from JCj{P) all 
unbounded faces, that is, faces that contain a vertex of F^o- 

If the defining inequalities of P are given, that is, P = {x : Ax < b}, the 
face lattice of P can be computed by an algorithm of Fukuda et al. |12) in 
time 0{rn-i{A, b)-ip), where i{A, b) is the time to solve a linear program with 
input size equal to the size of A and b ; its space complexity is 0{(p ■ logm + 
p{A,b)), where p{A,b) is the space needed to solve a linear program of the 
size of A and b. This algorithm outputs the faces as the sets of facets they 
are contained in. It is also possible to apply the algorithm to the unbounded 
polyhedron P. Note that the algorithm does not produce the Hasse diagram, 
but it can be computed using the algorithm in |20) . 

After the generation of one computes the set of vertices of F^o 

and removes the unbounded faces. Additional work is necessary if the faces 



COMPUTING BOUNDED SUBCOMPLEXES 



9 



of B{P) should be given by their vertex sets. This approach via the en- 
tire face lattice of P is not efficient when the bounded subcomplex B{P) is 
much smaller. It is not obvious to the authors whether or not the algorithm 
from |12) can be modified to compute the bounded faces only. 

The above algorithm avoids the explicit computation of the vertex-facet 
incidences (with unknown complexity) and leads to a polynomial total time 
algorithm if ip' ~ ip. In practice or if ip' <^ ip, Algorithm [1] might be faster. 

3.4. Polyhedra Given in Terms of Vertices and Rays. If the un- 
bounded polyhedron is given by the list of its vertices and rays, one can 
proceed analogously to the previous section by applying the algorithm of 
Fukuda et al. |12| to the dual of P. It is also easy to adapt their algorithm 
to work with vertices as input (their "restricted face of polyhedron" problem 
can also be solved via a linear program in this case). The faces are then 
given by their vertices and rays. This approach has the same drawbacks as 
the one in the previous section and can also be applied directly to P. 

The affine hull of each bounded k-iace of P is spanned hy k + 1 affinely 
independent vertices. For each (k + l)-tuple of the n vertices of P, one 
can solve one linear program to decide whether there exists a supporting 
hyperplane containing those k + 1 vertices and whose intersection with P is 
A;- dimensional. This immediately yields the following. 

Proposition 9. Let P be given in terms of vertices and rays, and let 5 be 
a fixed constant. Then the set of bounded faces up to dimension 5 can be 
computed in 0{n^ ■ £(V, const)) time. 

This algorithm does not produce the Hasse diagram of the (5-skeleton 
of B{P) directly, but this can be achieved via the algorithm in |20j . Al- 
ternative algorithms would result from computing the facets and the corre- 
sponding incidences and then applying the previously mentioned algorithms. 

Note that Algorithms [1] and [2] can easily be modified to produce the 5- 
skeleton as well. 

4. Simple Polyhedra 

In this section, we deal with the special case of simple polyhedra. A pointed 
d-dimensional polyhedron P is simple if each vertex is contained in precisely d 
facets. Note that P may not be simple, even if P is. Nevertheless, using 
a suitable generic construction one can guarantee that the corresponding 
polytope is simple, see [T^l Prop. 2.2]. 

For simple polyhedra. Algorithm [T] can be implemented more efficiently. 
Step E] can be performed in time 0{d ■ a), and the for- loop is executed at 
most 0{d ■ if') times, see |20) for more details. Note that the number of 
vertex-facet incidences is a = d -n. Thus, we obtain the following. 

Proposition 10. Given a simple polyhedron P and the vertex-facet inci- 
dences of P, the (Hasse diagram of the) bounded subcomplex of P can be 
computed in time 0{d • a • ip') = 0{d^ ■ n ■ ip') time. 

If a simple polyhedron P is given in terms of inequalities the reverse search 
algorithm by Avis and Fukuda [3] generates the vertices of P in 0{d-m) time 
per vertex, where d is the dimension of P. This is possible by performing 
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Input: inequality description of a simple, pointed, unbounded 

polyhedron P, far face F^q 
Output: (/o, ...,fd) = /-vector of B{P) 

1 compute the vertices and the vertex-edge graph F of P 

2 find c such that linear program maxjc^'^x : x € P} takes its maximum 
on and which is generic on vertices of P 

3 direct the edges of T along increasing c 

4 d dimP 

5 for k ^ d,d — 1, ... ,0 do 

6 hj. ^ number of vertices of P with out-degree k 

7 ^ number of vertices of Pqo with in-degree k 

8 lfk^T.Ui:u){h^-hT) 

Algorithm 3: Face numbers of simple polyhedra 



a ratio test that decides whether an edge is unbounded or not; see Avis [T] 
for details. During the vertex enumeration, the vertex-facet incidences of P 
and P can be computed on-the-fly at no extra cost. Using the incidences 
for P, we get the following. 

Corollary 11. Given a simple polyhedron P in inequality form, (the Hasse 
diagram of) B{P) can he computed in time 0(max{m, d ■ f'} ■ d-n). 

Kalai gave an algorithm to compute the vertex-facet incidences of a simple 
polytope from its vertex-edge graph |21| . This method can be modified 
to compute the f -vector (/o, /i, . . . , fd-i) of the bounded subcomplex of a 
simple d-polyhedron, where is the /c-th face number, that is, the number 
of faces of dimension k. Algorithm [3] presents this approach. 

Theorem 12. Algorithmic computes the face numbers of a simple polyhedron 
in 0{d -rn -n) time. 

Proof. The reverse search algorithm [3] produces the graph F of P in time 
0{d -rn-n). The running time of Step [1] dominates the remaining steps. 

Knowing the vertices of P and the far face is equivalent to knowing the 
vertices and rays of P. Each ray describes a direction in which P is un- 
bounded. Such a ray can be perturbed such that the corresponding linear 
objective function c takes distinct values on distinct vertices (see Step [2]). 

The numbers (/iq, hi, . . . , h^) computed in Step[6]form the h-vector of P, 
compare |21) . Its relationship with the /-vector, expressed in Step [HI is as 
follows: Each /c-face of P has a unique minimal vertex with respect to c. 
Conversely, each k-set of arcs which leave a fixed vertex v spans a fc-face 
such that V the minimum with respect to c is attained at v. We have to 
ignore the unbounded fc-faces of P, which are precisely the /c-faces of P 
whose maximum with respect to c is attained at some vertex of Pqo- D 

The formula in jl5( Prop. 2.4] for the relationship between the /- and 
/i-vectors of an unbounded polyhedron is wrong. The correct version is in 
Step [8] of Algorithm [31 
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5. Computational Results 



The following experiments were performed with the polymake system |13| . 
version 2.9.8. The hardware used was an AMD Athlon 64 X2 Dual Core 
Processor 4200 (4435.84 bogomips) with 4GB main memory running Debian 
Linux. We tested Algorithm [1] requiring the vertex-facet incidences of P as 
input only. The timings for the convex hull computations required are not 
given since discussing the various choices is a topic of its own; see [21 117| . 
Here we are focusing on the combinatorial aspects. Usually the convex hull 
computation takes much less time than the computation of the face lattice. 

The goals of our computations are threefold. 

o The limits of our algorithm can be estimated. 

o Several examples for which the computation of B{P) is interesting are 

investigated and corresponding results are presented, 
o A rough estimation of the sizes up to which a computation of B{P) is 

sensible, independent of the approach, can be derived (for our examples). 

Our theoretical analysis so far was based on the assumption that sets 
are stored as bitsets; polymake offers a suitable data type (wrapping an 
implementation of the CMP |14j). However, the representation of sets via 
balanced trees seems to be superior in a typical scenario. This is the one 
used in the tests below. 

5.1. Dwarfed Cubes. The dwarfed d-cuhe is the polytope 



The polytope D has m = 2d+ 1 facets and n = + 1 vertices. Moreover, 
a = d -n as the polytope is simple. The interest in these polytopes comes 
from the fact that the dwarfed cubes provide difficult input to some classes 
of convex hull algorithms, see Avis et al. [2]. To produce an unbounded 
polyhedron D, we send the dwarfing facet ^ Xj = 3/2 to infinity by reversing 
the construction from Proposition [TJ The dwarfing facet contains d{d — 1) 
vertices; hence the bounded subcomplex has n = + 1 — d{d — 1) = d + 1 
vertices and f3 = n- d = d'^ + d incidences between vertices and facets. The 
bounded subcomplex B{D) is a star-like graph with d+1 nodes (and d edges); 
that is, the number (p' of bounded faces equals 2d + 2, including the empty 
face. 

Table [1] contains the results of the computations. They show that B{P) 
could be computed up to a very high dimension (> 70) with a moderate 
computing time; the key reason for this seems to be the small number of 
bounded faces (p' . 

Remark 13. In order to count the faces of the (unbounded) dwarfed cube D 
we can employ Algorithm [3l To this end consider the linear objective function 
^ Xi which is generic on the unbounded dwarfed cube, and which takes 
its maximum on the dwarfing facet. This linear objective function gives a 
direction on each bounded edge or ray of D. The origin is the unique node of 
out-degree d, making up for 2"^ non-empty faces of D whose minimum with 
respect to ^ Xi is 0. The d neighbors of 0, the d unit vectors, are nodes of 
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Table 1. Results for Dwarfed cubes 



d 


771 


71 






timp ( SI 


5 


11 


26 


130 


12 


0.07 


10 


21 


101 


1010 


22 


0.11 


15 


31 


226 


3390 


32 


0.42 




41 

J. 






49 


J- .O 1 


25 


51 


626 


15650 


52 


6.89 


30 


61 


901 


27030 


62 


28.35 


35 


71 


1226 


42910 


72 


72.93 


40 


81 


1601 


64040 


82 


159.64 


45 


91 


2026 


91170 


92 


324.86 


50 


101 


2501 


125050 


102 


593.73 


55 


111 


3026 


166430 


112 


1039.30 


60 


121 


3601 


216060 


122 


1743.58 


65 


131 


4226 


274690 


132 


2811.10 


70 


141 


4901 


343070 


142 


4457.96 


75 


151 


5626 


421950 


152 


6823.86 



out-degree d—1, making up for an additional (i-2'^ ^ non-empty faces. Thus 
ip{D) = 2'^ + d- 2'^-^ + 1 (including and D). 

The poset J'(D) contains 2"^ — 1 elements for faces of the first kind above 
(not containing the face D itself) and d elements corresponding to the d unit 
vectors. Thus ^p" = 2'^ -\- d (including 0). In contrast, ip' is only linear in d. 

5.2. Tight Spans of Metric Spaces. The tight span (or injective hull) Tm 
of a finite metric space M : [d\ x [d] — )• R (see Dress |11| and Isbell |16) ) is 
defined as the bounded complex of the polyhedron 

Pj^l = |x G E'' : Xj > M{i,j) for all 1 < j < d| . 

If M is generic enough, as in our examples, all inequalities define facets, hence 
_ _ d{d+i) _|_ ^ (including i^oo)- It was remarked by Sturmfels and Yu |25| 
that Tm is dual to the complex of inner faces of the regular subdivision of 
the second hypersimplex 

A(2, d) := conv {ej + Cj : 1 < « < i < d} , 

obtained from interpreting M as a height function. Tight spans are relevant 
for applications in algorithmic biology, more precisely in phylogenetics, cf. 
Dress et al. 

Two ways to obtain special metric spaces and corresponding examples are 
described next. 

5.2.1. Thrackle Metric. A special triangulation of A(2, d), called the thrackle 
triangulation was introduced by Stanley |24j and thoroughly investigated by 
De Loera et al. [8]. It turns out that the corresponding metric maximizes 
the number of faces of the tight span |15| and is equivalent to the tight span 
of the maximal circular split system, see Bandelt and Dress [H Section 3]. 
From |15| Theorem 5.5], we deduce that n = 2'^~^ -\- d and 

= 1 ((1 - V2)'' + (1 + ^/2)'^) + 1. 
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Table 2. Results for thrackle metrics 





d 


m 


n 


a 




time (s) 




3 


7 


7 


24 


8 


0.01 




4 


11 


12 


60 


18 


0.01 







16 


21 


135 


42 


n m 

U.Ui 




6 


22 


38 


288 


100 


0.03 




7 


29 


71 


602 


240 


0.15 




8 


37 


136 


1248 


578 


0.96 




9 


46 


265 


2565 


1394 


7.38 




10 


56 


522 


5210 


3364 


61.90 




11 


67 


1035 


10450 


8120 


559.08 




12 


79 


2060 


20736 


19602 


5239.04 




13 


92 


4109 


40820 


47322 


54302.46 






Table 3. Results for random metrics 




d 


m 


n 


a 




time (s) 


stddev (s) 


5 


16 


21.00 


135.00 


42.00 


0.03 


0.00 


6 


22 


37.99 


287.94 


99.92 


0.06 


0.00 


7 


29 


70.65 


599.55 


237.20 


0.19 


0.01 


8 


37 


134.95 


1247.60 


568.96 


1.17 


0.09 


9 


46 


261.94 


2609.46 


1365.28 


9.52 


0.99 


10 


56 


513.55 


5495.50 


3275.68 


86.69 


9.69 


11 


67 


1008.46 


11588.06 


7802.88 


841.02 


109.44 


12 


79 


1997.28 


24627.36 


18709.52 


9043.11 


1351.32 



Note that the right hand side is always integral. 

Table [2] shows the results. Not surprisingly, it turns out that with increas- 
ing dimension the the computation time drastically increases along with the 
number of bounded faces. 

5.2.2. Random Metrics. Our random metrics are obtained by taking the dis- 
tances M(i,j) to be uniformly distributed in the interval [1,2]. Since this is 
generic (with probability 1), we have m = '^^'^^^^ + 1. The sample size for 
each dimension where 100 metrics. For n, a, and ip' , we state the arithmetic 
mean, and for the computation time the mean together with the standard 
deviation. 

Table [3] presents the results. Compared to the thrackle metrics (Table [2|), 
random metrics have fewer bounded faces, but their computation times are 
larger. 

5.3. Tropical Polytopes. Let V = [vik) be an s x t-matrix with real coef- 
ficients. We define the polyhedron 

Ev ■■= {{u,w) : Ui+Wk<Vik] , 

where w E and w E M*. Considering {u,w) with sufficiently small coordi- 
nates, one can see that Ey is not empty. Moreover, if (u, w) E Ey then for 
all A E M we have {u + Xl,w — XI) = {u,w) + X{l,—l) E Ey- Hence the one- 
dimensional subspace M(l, —1) is contained in the lineality space of Ey, and 
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Table 4. Results for tropical cyclic polytopes 



(s,t) 


d 


m 


n 


a 


^' 

r 


time (s) 


(3,3) 


5 


10 


12 


72 


14 


0.04 


/a a\ 

(4,4) 


7 


17 


28 


244 


64 


0.04 


(5,5j 


n 

y 


ZD 


oU 


o4U 


OZZ 


n /in 
U.4U 


(6,6) 


11 


37 


264 


3144 


1684 


17.52 


(7,7) 
(8,8) 


13 


50 


938 


12614 


8990 


1198.35 


15 


65 


3448 


52392 


48640 


139091.23 


(3,10) 


12 


31 


68 


1003 


182 


0.22 


(3,20) 


22 


61 


233 


5903 


762 


9.17 


(3,30) 


32 


91 


498 


17703 


1742 


110.16 


(3,40) 


42 


121 


863 


39403 


3122 


814.14 


(3,50) 


52 


151 


1328 


74003 


4902 


4418.48 


(3,60) 


62 


181 


1893 


124503 


7082 


15595.14 


(3,65) 


67 


196 


2213 


156653 


8322 


26858.18 


(3,70) 


72 


211 


2558 


193903 


9662 


44651.71 



hence we can consider Ey as a polyhedron in the quotient ]R*"'"*/IR(1, — 1). 
The polyhedron Ey in M*"'"*/M(l, —1) is pointed, and projecting its bounded 
subcomplex to ffi* (or, alternatively, to M") yields the tropical polytope de- 
fined by V; see Develin and Sturmfels 0. The bounded subcomplex of Ey 
is dual to the regular subdivision of the product of simplices As-i x Af_i, 
obtained by interpreting the matrix F as a lifting function. Here A^. denotes 
a simplex of dimension r. 

We now translate the parameters for tropical polytopes into the parame- 
ters that we used in our algorithms. The dimension of the polyhedron Ey 
(in the quotient) equals d = s + 1 — 1. Its number m of facets is less than or 
equal to s ■ t. Throughout we have m = m + 1. The number n of vertices 
satisfies n < ('^^1^^). 

5.3.1. Tropical Cyclic Polytopes. The tropical cyclic polytope with parame- 
ters (s, t) is given by the matrix V = {vik) with = i-k. The corresponding 
subdivision of As_i x At-i is known as the staircase triangulation; see Block 
and Yu [6]. The polyhedron Ey is simple in this case. 

5.3.2. Tropical Permutohedra. Each permutation a on the t numbers from 
tot— \ can be identified with the vector ((t(0), (t(1), . . . , (T{t — 1)). This way 
each permutation contributes one row of a (t!) x t-matrix V . We call the cor- 
responding tropical polytope a tropical permutohedron. The polyhedron Ey 
is not simple for t > 3. 



Table 5. Tropical permutohedra 



(s,t) 


d 


m 


n 


a 




time (s) 


(6,3) 


8 


19 


24 


261 


50 


0.05 


(24,4) 


17 


97 


152 


6532 


1424 


9.07 


(120,5) 


124 


601 


1420 


276725 


76282 


143535.58 
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Table H] shows the results for tropical cyclic polytopes and Table [5] for 
tropical permutohedra. Both cases show a steep increase in computation 
time with increasing dimension. The remarkable fact is that examples of 
these sizes can be handled at all. 

6. Concluding Remarks and Open Questions 

We presented combinatorial algorithms to compute (the Hasse diagram of) 
the bounded faces of an unbounded pointed polyhedron. Algorithm [1] which 
takes the vertex-facet incidences of P as input, was also shown to work 
in practice via extensive computations. It seems that the examples which 
we presented in the last section show the limits of what one can currently 
compute in acceptable time, unless some (possibly fundamentally) different 
idea comes up. For instance, a central open question is the following. 

Question. Is there a polynomial total time algorithm to compute B{P) for 
general P, given the inequalities? 

For instance, this would follow if it were possible to modify the algorithm 
of Fukuda et al. |12) to compute the bounded faces only. 

In order to get an idea about the size of a bounded subcomplex it would 
be interesting if it were possible to compute statistical information without 
generating all faces. 

Question. Is there a polynomial time algorithm to compute the f -vector of 
L{P) orB{P)? 

See the discussion in Section [J] for simple polyhedra. 
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